DESY 97-183 



Observation of Scaling Violations 
in Scaled Momentum Distributions 

at HERA 

ZEUS Collaboration 



Abstract 

Charged particle production has been measured in deep inelastic scattering (DIS) 
events over a large range of x and using the ZEUS detector. The evolution of the 
scaled momentum, Xp, with Q^, in the range 10 to 1280 GeV^, has been investigated in 
the current fragmentation region of the Breit frame. The results show clear evidence, in 
a single experiment, for scaling violations in scaled momenta as a function of Q^. 
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1 Introduction 



The observation of scaling violations in structure functions measured in neutral current, 
deep inelastic scattering (DIS) helped to establish Quantum Chromodynamics (QCD) as the 
theory of strong interactions and has led to measurements of the strong coupling constant, 
as- Similar scaling violations are predicted in the fragmentation functions, which represent the 
probability for a parton to fragment into a particular hadron carrying a fraction of the parton's 
energy. Fragmentation functions incorporate the long distance, non-perturbative physics of the 
hadronization process in which the observed hadrons are formed from final state partons of the 
hard scattering process and, like structure functions, cannot be calculated in perturbative QCD, 
but can be evolved from a starting distribution at a defined energy scale. If the fragmentation 
functions are combined with the cross sections for the inclusive production of each parton 
type in the given physical process, predictions can be made for scaling violations in the scaled 
momentum spectra of final state hadrons [0]. These scaling violations allow a measurement of 
as and such studies have been performed at LEP P, Q by incorporating lower energy PETRA 
data. 

The event kinematics of DIS are determined by the negative square of the 4-momentum transfer 
at the positron vertex, = — g^, and the Bjorken scaling variable, x = Q^/2P ■ q, where P 
is the four-momentum of the proton. In the Quark Parton Model (QPM), the interacting 
quark from the proton carries the four-momentum xP. The variable y, the fractional energy 
transfer to the proton in its rest frame, is related to x and by ?/ ~ Q'^/xs, where ^/s is 
the positron-proton centre of mass energy. A natural frame in which to study the dynamics 
of the hadronic final state in DIS is the Breit frame H, which has been used in previous 
studies of QCD effects ||^, ||, |^ at HERA. In this frame the exchanged virtual boson is purely 
space-like with 3-momentum q = (0,0, —Q), the incident quark carries momentum Q/2 in the 
positive Z-direction, and the outgoing struck quark carries Q/2 in the negative Z-direction. 
A final state particle has momentum in this frame, and is assigned to the current region 
if Pz is negative, and to the target frame if is positive. The advantage of this frame lies 
in the maximal separation of the outgoing parton from radiation associated with the incoming 
parton and the proton remnant, thus providing the optimal environment for the study of its 
fragmentation. 

In this analysis the inclusive charged particle distributions of the scaled momenta, Xp, in the 
current region of the Breit frame are measured; Xp is the momentum of a track measured in 
the Breit frame, p^, scaled by Q/2, the maximum possible momentum (ignoring effects due to 
the intrinsic of the quark within the proton). 

The scaled momentum distributions are studied as a function of Q and x, in the range 6 x 10~^ < 
X < 5 X 10~^ and 10 < < 1280 GeV^. Thus the evolution of fragmentation functions can 
be observed within a single experiment over a wide range of Q- A similar analysis has recently 



been performed by the HI collaboration [jTO|. 
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2 Experimental Setup 



The data presented here were taken in 1994 at the positron-proton colhder HERA using the 
ZEUS detector. During this period HERA operated with positrons of energy Ee = 27.5 GeV 
and protons with energy 820 GeV. The data of this analysis correspond to a luminosity of 



2.50 ± 0.04 pb . A detailed description of the ZEUS detector can be found in |11, O; we 



present here a brief description of the components most relevant to the present analysis. 

Throughout this paper we use the standard ZEUS right-handed coordinate system, in which 
X = Y = Z = Ois the nominal interaction point, the positive Z-axis points in the direction of 
the proton beam (referred to as the forward direction) and the X-axis is horizontal, pointing 
towards the centre of HERA. 

The tracking system consists of a vertex detector (VXD) [jl3[ and a central tracking cham- 
ber (CTD) enclosed in a 1.43 T solenoidal magnetic field. Immediately surrounding the 
beampipe is the VXD which consists of 120 radial cells, each with 12 sense wires. The CTD, 
which encloses the VXD, is a drift chamber consisting of 72 cylindrical layers, arranged in 9 
superlayers. Superlayers with wires parallel to the beam axis alternate with those inclined at a 
small angle to give a stereo view. The single hit efficiency of the CTD is greater than 95% and 
the measured resolution in transverse momentum for tracks with hits in all the superlayers is 
ap^/pT = 0.005pT 00.016 {pr in GeV). 



Outside the solenoid is the uranium-scintillator calorimeter (CAL) |jT5[, which is divided into 
three parts: forward, barrel and rear covering the polar regions 2.6° to 36.7°, 36.7° to 129.1° and 
129.1° to 176.2°, respectively. The CAL covers 99.7% of the solid angle, with holes of 20 x 20 
cm^ in the centres of the forward and rear calorimeters to accommodate the HERA beam pipe. 
Each of the calorimeter parts is subdivided into towers which are segmented longitudinally into 
electromagnetic (EMC) and hadronic (HAG) sections. These sections are further subdivided 
into cells each of which is read out by two photomultipliers. From test beam data, energy 
resolutions of ge/E = 0.18/v^ for electrons and (Je/E = 0.35/V^ for hadrons {E in GeV) 
have been obtained. 

The small angle rear tracking detector (SRTD) [|l6l, which is attached to the front face of the 
rear calorimeter, measures the impact point of charged particles at small angles with respect 
to the positron beam direction. 

The luminosity is measured by means of the Bethe-Heitler process ep — >■ e7p, by detecting the 
photon in a calorimeter [|T^] positioned at Z = —107 m which has an energy resolution of aE/E 

= 0.18/ E{GeV) under test beam conditions. The luminosity calorimeter is also used to tag 
photons from initial state radiation in DIS events. 



3 Kinematic Reconstruction 



The ZEUS detector is almost hermetic, allowing the kinematic variables x and to be recon- 
structed in a variety of ways using combinations of positron and hadronic system energies and 
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angles. Variables calculated only from the measurements of the energy, E'^, and angle, 6^, of the 
scattered positrons are denoted with the subscript e, whilst those calculated from the hadronic 
system measurements, with the Jacquet Blondel method |T^, are denoted by the subscript JB. 
Variables calculated by these approaches are only used in the event selection. In the double 
angle method [0, denoted by DA, the kinematic variables are determined using 6e and the 
angle (which is the direction of the struck quark in QPM), defined from the hadronic final 
state: 



{J:^Pxy + {J:^PYr + {J:^{E-pz)y ' 

where the sums run over all CAL cells i, excluding those assigned to the scattered positron, 
and p = {px,Py,Pz) is the 3-momentum assigned to a cell of energy E. The cell angles are 
calculated from the geometric centre of the cell and the vertex position of the event. This angle 
is then combined with the measured angle of the scattered positron to calculate x and Q^. The 
two angles can also be used to recalculate the energy of the scattered positron: 



^, ^ 2^esin7„ 

sin 7^ + sin 6e — sin(7^ + 6^) 

The DA method was used throughout this analysis for the calculation of the boosts and the 
kinematic variables because it is less sensitive to systematic uncertainties in the energy mea- 
surement than the other methods discussed above. 

An additional method of measuring y and (the PT method |^) was used as a systematic 
check of the kinematic reconstruction to determine the binning of the data and the boost vector. 
This method uses a more sophisticated combination of the information from the measurements 
of both the hadronic system and the positron. 

In order to boost to the Breit frame it is necessary to calculate the velocity of the Breit frame 
with respect to the laboratory frame, which is given by /? = {q + 2xP) / (go + 2xPo) where 
(go, <?) and {Pq, P) are the 4-momenta of the exchanged photon and the incident proton beam 
respectively. The event is then rotated so that the virtual photon is along the negative Z- 
axis. The Breit frame boost vector was reconstructed using -E^^ and the polar and azimuthal 
angles measured from the impact point of the scattered on the calorimeter or SRTD. The 
4-momentum vectors of the charged particles were boosted to the Breit frame, assuming the 
pion mass to determine the particle's energy, and were assigned to the current region if pf < 0. 



4 Event Reconstruction and Selection 



The triggering and online event selections were identical to those used for the measurement of 
the structure function F2 
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To reduce the effects of noise due to the uranium radioactivity on the calorimetric measurements 
all EMC (HAC) cells with an energy deposit of less than 60 (110) MeV were discarded from 
the analysis. For cells with isolated energy deposits this cut was increased to 100 (150) MeV. 
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The track finding algorithm starts with hits in the outer axial superlayers of the CTD. As the 
trajectory is followed inwards towards the beam axis, hits from inner superlayers are added to 
the track. The momentum vector is determined in a 5-parameter helix fit. The reconstructed 
tracks used in this analysis are associated with the primary event vertex and have pt > 150 MeV 
and \ri\ < 1.75, where 77 is the pseudorapidity given by — ln(tan(^/2)) with 6 being the polar 
angle of the measured track with respect to the proton direction. This is a region of high CTD 
acceptance where the detector response and systematics are best understood. 

Further selection criteria were applied both to ensure accurate reconstruction of the kinematic 
variables and to increase the purity of the sample by eliminating background from photopro- 
duction processes. These cuts were: 

• E'^> 10 GeV, to achieve a high purity sample of DIS events; 

• QIa > 10 GeV2; 

• 2/e < 0.95, to reduce the photoproduction background; 

• yjB > 0.04, to give sufficient accuracy for DA reconstruction; 

• 35 < 6 = J2i {E — pz) < 60 GeV summed over all calorimeter cells, to remove photopro- 
duction events and events with large radiative corrections. 

• |X| > 16 cm or |y| > 16 cm, where X and Y are the impact position of the positron on 
the GAL as determined using the SRTD to avoid the region directly adjacent to the rear 
beam pipe. 

• —40 < ^vertex < 50 cm, to reduce background events from non ep collisions. 

In total, 68066 events satisfy the above cuts and are reconstructed in the {x, Q'^) bins (as 
calculated by the DA method) that are listed in Table |I]. The sizes of the bins in x and 
were chosen to have good statistics in each bin and to limit the migrations between bins 0]. 



5 QCD Models and Event Simulation 



Monte Carlo event simulation is used to correct for acceptance and resolution effects. The detec- 



tor simulation is based on the GEANT 3.13 |21| program and incorporates our best knowledge 
of the apparatus. 

To calculate the acceptance, neutral current DIS events were generated, via the DJANGO 
program [Q, using HERACLES |^ which incorporates first order electroweak corrections. 
The QCD cascade was modelled with the colour-dipole model including the boson-gluon fusion 
process, using the ARIADNE 4.08 program. In this model coherence effects are implicitly 
included in the formalism of the parton cascade. This program uses the Lund string fragmen- 
tation model 1^ for the hadronisation phase, as implemented in JETSET 7.3 PBI. Two Monte 



Carlo samples were generated, 1.6 pb"^ with > 3 GeV^ and 2.6 pb"^ with > g QeV^, 
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using the MRSA parameterisation of the parton distribution functions. Another approach 
to modelhng the parton cascade is included in the LEPTO(6.5_l) |^ program, which incor- 
porates the LO matrix element matched to parton showers (MEPS). This program also uses 
the Lund string fragmentation model and was used for the generator level calculations used in 
comparisons to our data. 

For the studies of the systematics, two additional samples of events were generated (1.2 pb~^ 
with > 4 GeV^ and 3.2 pb"^ with > iqo GeV^) using the HERWIG 5.8c Monte Carlo |29|, 
where no electroweak radiative corrections were applied. In HERWIG, coherence effects in the 
QCD cascades are included by angular ordering of successive parton emissions and a clustering 
model is used for the hadronisation |30, 31]. For this sample the parameterisation of the parton 
distribution functions was the MRSA' set [0 . The MRSA and MRSA' parameterisations have 
both been shown to describe reasonably well the HERA measurements of the proton structure 
function F2 in the (x, Q"^) range of this analysis [p3, . 



6 Correction Procedure 



Monte Carlo event generator studies were used to determine the mean charged particle accep- 
tance in the current region as a function of (x, Q^). The chosen analysis intervals in (x, Q^) 
correspond to regions of high acceptance, 74 to 96%, in the current region of the Breit frame. 



xda range 


Ql)A (GeV^) range 


No. of events 


6.0 - 12.0 X 10"^ 


10-20 


13898 


1.2 - 2.4 X 10-3 


10-20 


11899 




20-40 


8484 




40-80 


5093 


2.4 - 10.0 X 10-3 


20-40 


9399 




40-80 


9493 




80 - 160 


5031 




160 - 320 


1369 


1.0 - 5.0 X 10"^ 


160 - 320 


2131 




320 - 640 


916 




640 - 1280 


353 



Table 1: The {x,Q'^) analysis bins showing the the numbers of accepted events in the bins as 

reconstructed by the DA method. 



Uncertainty in the reconstruction of the boost vector, (3, was found to be the most signif- 
icant factor on the resolution of Xp and it leads to the choice of variable bin width in Xp. 
Migration of tracks from the current region to the target region was typically ~ 6 to 8% of 
the tracks generated in the current region. Migrations into the current region from the target 
fragmentation region are typically less than 5% of the tracks assigned to the current region for 
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> 320 GeV^. For 10 < < 320 GeV^ these migrations are generally of the order of 10 to 
15% . At < 40 GeV^ and low values of y this assignment can be as high as ~ 25%, since in 
the low y region the hadronic activity is low and the measurement of becomes distorted by 
noise in the calorimeter leading to a worse x resolution and hence an uncertainty in (3. 

The correction procedure is based on the detailed Monte Carlo simulation of the ZEUS detector 
with the event generators described in the previous section. Since the ARIADNE model gives 



the best overall description of our observed energy flow it is used for the standard corrections 
to the distributions. 

The data are corrected for trigger and event selection cuts; event migration between (x, Q"^) 
intervals; QED radiative effects; track reconstruction efficiency; track selection cuts in and 
77; track migration between the current and target regions; and for the products of A and 
decays which are assigned to the primary vertex. 

Correction factors were obtained from the Monte Carlo simulation by comparing the generated 
distributions, without A and Kg decay products, with the reconstructed distributions after the 
detector and trigger simulations followed by the same reconstruction, selection and analysis as 
the real data. The correction factors, F{x^, were calculated for each Xp bin with the formula: 




P/ obs 



where Agen (Aobs) is the number of generated (observed) Monte Carlo events in each (x, Q"^) 
interval and n is the number of charged particles (tracks) in the current region in the corre- 
sponding (x, Q^) interval. The correction factors are typically between 1.0 and 1.3, except in 
the lowest Xp bin (0.02 < Xp < 0.05) of the two lowest bins where they are about 1.5. 



7 Results 

The inclusive charged particle distributions, 1/ atotda / dxp where atot is the DIS cross section 
in the chosen (x, Q^) bin, are shown in figure |l|. In the low (x, Q'^) regions these distributions 
peak at Xp ~ 0.2 with a broad tail out to 1. As increases the data fall off more sharply 
from the minimum Xp. In the same figure the ARIADNE Monte Carlo generator predictions 
are also shown, which are in good agreement with the data. The increasingly steep fall-off 
towards higher values of Xp as increases corresponds to the production of more particles 
with a smaller fractional momentum, and is indicative of scaling violation in the fragmentation 
function. These scaling violations can be seen more clearly if the data are plotted in bins of 
fixed Xp as a function of Q^. Figure |^ shows the charged particle distributions with statistical 
errors combined in quadrature with the systematic errors which will be discussed below. For 
> 80 GeV^ the distributions rise with at low Xp and fall-off at high Xp and high Q^. By 
measuring the amount of scaling violation one can ultimately measure the amount of parton 
radiation and thus determine as- Below = 80 GeV^ the fall-off is due to depopulation of the 
current region discussed later. 
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7.1 Systematic Studies 



The systematic uncertainties in the measurement can be divided into three classes: errors due 
to event reconstruction and selection; to track reconstruction and selection; and to the Monte 
Carlo generator used. 

• Event Reconstruction and Selection 

The event selection procedure was checked by altering the ?/e, Vjb and S cuts: the resulting 
shifts of the corrected distributions were small, with all points moving systematically in 
the same direction at the ±1 to 3% level. By removing the noise suppression, described in 
section ^, a systematic shift at the ±1 to 3% level is produced. Using the boost calculated 
from the PT method produced a larger shift: for < 160 GeV^ the systematic shift is 
+ 10 to 15% at large Xp, whilst +2 to 5% at small x^; for > 160 GeV^ the systematic 
shift is constant with Xp at the +2 to 5% level. 

• Track Reconstruction and Selection 

The major systematic on the track reconstruction was obtained considering all recon- 
structed tracks as opposed to only those assigned to the primary vertex. This produced 
a typical shift of 2 to 6% with no systematic dependence on the value of Xp, except in 
the bin 0.02 < Xp < 0.05 at low and x where the systematic shift was up to +15%. 
Tightening the tracking cuts on \ri\ and pt particularly affected those (x, Q^) bins where 
the track acceptance is strongly dependent on the negative t] cut. The general trend is for 
the cross section to be shifted to higher values. The bins most affected have an average 
shift of up to 5%, whilst in general the shift is <,1%. 

• Monte Carlo Generator 

Using a different Monte Carlo generator (HERWIG rather than ARIADNE) led to dis- 
tributions which were systematically lower by about 10% in the range of < Xp < 0.3. 
In the range Xp > 0.5, the HERWIG corrected distributions were systematically higher 
(lower) by 5 to 10 % for < 40 GeV^ {Q^ > 100 GeV^). In the range 0.3 < Xp < 0.4 
the two generators gave results that were in good agreement. 

All positive (negative) systematic shifts in each of the Xp bins were combined in quadrature to 
give an estimate of the overall positive (negative) systematic uncertainty on the measurement. 
These systematic shifts do not affect the observation of scaling violations. 



7.2 Discussion 



Fragmentation in DIS of a quark carrying momentum Q/2 in the Breit frame may be compared 
to fragmentation in e~^e~ annihilation of the produced q and q, each carrying momentum 
y/Se+e-/2. In figure the ZEUS data are compared at = Sg+e- to e+e~ data [^, divided by 
two to account for the production of both a q and q. The e+e" data have also been corrected 
by ~ 8% for the decay products of A and Kg, using the JETSET 7.3 Monte Carlo tuned to 
DELPHI data. In the range shown there is good agreement between the current region of 
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the Breit frame in DIS and the e^e experiments. When our data are compared with the lower 
energy SPEAR |^ data (i/i = 5.2, 6.5 GeV, not shown) discrepancies begin to show up. They 



can be explained in terms of the kinematic depopulation of the current region described in the 
next paragraph. 

The turnover observed in the ZEUS data at low Xp and low Q^, figure ||, can be attributed 
to processes not present in e"'"e~ (eg scattering off a sea quark and/or boson gluon fusion) 
depopulating the current region p. This is best illustrated when discussing the production 
of a pair of partons in DIS with a large invariant mass, s When ^ s, the radiation 
is emitted in the direction of the struck quark in the QPM. However, at low and low x, s 
is likely to be bigger than and the radiation will be emitted in the direction opposite to 
the QPM direction and will therefore not be in the current region as defined in section |l|. In 
particular, the boson-gluon fusion diagram, which dominates at low-x, provides a significant 



cross section for large mass radiation |^, thus producing this depopulation. 



In figure |^ the data are compared to two leading-log Monte Carlo predictions which are im- 
plemented within the Lund fragmentation framework. Our data are well described by the 
ARIADNE generator, colour dipole model (CDM), over the full range of Q^. This is not true 
for the LEPTO, matrix element +parton showers (MEPS), model. In particular LEPTO over- 
estimate the data at low Xp and there is a greater dependence than that suggested by the 
data. A feature of both Monte Carlo models is a trend that, at fixed , the charged particle 
distribution increases with x. Such a trend is also observed in the data, see fig |. No tuning of 
the Monte Carlo parameters has been performed to improve the agreement with the data. 

Our results can be compared to the next-to-leading order (NLO) QCD calculations, as imple- 
mented in CYCLOPS [^, of the charged particle inclusive distributions in the restricted region 
> 80 GeV^ and Xp > 0.1, where the theoretical uncertainties are small and unaffected by the 
hadron mass effects which are not included in the fragmentation function [40]. This comparison 
is shown in figure ^ The NLO calculation combines a full next-to-leading order matrix element 
with the MRSA' parton densities (with Aqcd = 230 MeV) and NLO fragmentation functions 



derived by Binnewies et al. from fits to e"*"e data [^. The data and the NLO calculations 
are in good agreement, supporting the idea of universality of quark fragmentation. 



8 Conclusions 

Charged particle distributions have been studied in the current region of the Breit frame in 
DIS over a wide range of values. These results show clear evidence in a single experiment 
for scaling violations in scaled momenta as a function of Q^. The data are also well described 
with NLO calculations. 

The comparison between e'^e" data at = Se+e- and the current region of the Breit frame 
in DIS for > 80 GeV^ shows also a good agreement. The observed charged particle spectra 
are consistent with the universality of quark fragmentation in e~^e~ and DIS. 
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ZEUS 1994 




X 6.0-12.0x10"* 1.2-2.4x10'^ 2.4-10.0x10'^ 1.0-5.0x10'^ 



Figure 1: The inclusive charged particle distributions, l/atot da/dxp^ in the current fragmentation 
region of the Breit frame compared with the generated distributions from the Monte Carlo (ARIADNE 
4.08). The statistical errors are generally within the size of the points. 
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Figure 2: The inclusive charged particle distribution, 1/atot da/dxp, in the current fragmentation 
region of the Breit frame. The filled data points are from ZEUS. The full errors are statistical and 
systematic combined in quadrature. The thick lines denote the statistical error. The open points 
represent data from e+e" experiments divided by two to account for q and q production (also corrected 
for contributions from Kg and A.) 
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Figure 3: The inclusive charged particle distribution, 1/crtot da/dxp, in the current fragmentation 
region of the Breit frame. The data points are from ZEUS. The full errors are statistical and systematic 
combined in quadrature. The thick lines denote the statistical error. The curves represent leading-log 
Monte Carlo models: the full line is ARIADNE (CDM) and the dashed line is LEPTO (MEPS). For 
clarity of presentation, only the higher x bin is shown in each interval. 
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Figure 4: The inclusive charged particle distribution, 1/atot da/dxp, in the current fragmentation 
region of the Breit frame compared to the NLO calculation, CYCLOPS . For clarity of presentation, 
only the higher x bin is shown in each interval. 



15 



